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system.  Ananev,  et  al.,  have  recently  developed  numerical 

3  ll 

analysis  on  the  efficiency  of  a  telescopic  resonator'”  . 
However,  no  tractable  results  are  available  for  the  achieve¬ 
ment  of  the  resonator  design  due  to  the  nonlinear  behavior 
of  the  rate  equation  and  the  coupled  optical  elements  of  the 
system. 

In  the  present  research  report,  we  will  analyze  the 
optimum  condition  for  the  energy  extraction  from  confocal, 
positive  branch,  unstable  resonators  for  Exclmer  lasers. 
Both  cylindrical  and  spherical  mirrors  for  20  and  30  gain 
elements,  respectively,  will  be  investigated.  For  tractable 
results,  analysis  on  simplified  systems  will  be  investigated 
and  compared  with  the  exact  numerical  results  for  more  real¬ 
istic  systems.  The  geometric-optics  approximation,  which 
has  been  shown  to  be  valid  for  systems  with  large  Fresnel 
numbers^*®,  will  be  employed,  i.e.,  the  small  deviation  of 
the  eigenvalue  due  to  the  multipass  diffraction  loss  and  the 
non-uniform  phase  profile  oh  the  output  mirror  (within  ten 
degrees  for  Fresnel  number  >  50)  will  be  Ignored. 

The  transverse  and  horizontal  variation  of  the 
forward  and  backward  Intensity  will  be  numerically  calculat¬ 
ed  for  the  case  of  constant  unsaturated  gain.  Analytic 
results  based  on  the  mean-field  approximation  will  be  deriv¬ 
ed  and  compared  with  the  numerical  results.  The  beam  quali¬ 
ty  based  on  the  focusability  of  the  far-field  will  be  ana¬ 
lyzed.  Finally,  some  suggested  measurements  which  provide 
system  parameters  information  such  as  optimal  magnification 
and  maximum  output  power  and  the  gain/loss  are  discussed. 

There  are  five  appendices  in  this  report  including 
the  computer  code. 


II 


THE  SYSTEM  AID  THE  PROBLEMS 


To  obtain  tha  energy  extraction  from  a  confocal 
unstable  resonator,  as  shown  in  Figure  1,  the  intensity 
inside  the  loaded  cavity  whose  gain  function  is  usually 
given  by  a  homogeneous  broadening  should  be  calculated. 
Subjected  to  the  oscillation  threshold  condition  (loss  - 
gain),  the  maximum  energy  extraction  may  be  obtained  by 
optimizing  system  parameters  such  as  reflectivity  and  magni¬ 
fication.  As  shown  by  the  flowchart  in  Figure  2,  our  final 
goal  is  to  find  the  resonator  efficiency  defined  by  [output 
power ]/Cgain] ,  where  the  output  power  is  given  by 

*  -  5  /  lout*8'-  C,] 

*  *  a/M  * 

Here  IQUt  is  the  output  intensity  evaluated  at  z  •  L  and  is 

integrated  over  the  outcoupling  regime  and  normalized  by  S; 

2 

dS*  -  2  irrdr,  S  -  va  for  spherical  mirrors  and  dS'  -  dr,  S 

I 

•  a  for  cylindrical  mirrors. 

In  general,  the  output  intensity  cannot  be  solved 
analytically  from  the  rate  equation  and  hence  computer  simu¬ 
lation  is  needed  for  the  resonator  efficiency.  We  note  that 

the  conventional  analysis  based  on  P  •  (1  -  M*n)  I  .  over- 

out 

simplifies  the  system  and  the  assumption  of  a  constant  IQut 

may  cause  considerable  error  in  the  prediction  of  the  energy 

>  extraction.  Before  showing  our  new  technique  for  obtaining 

resonator  efficiency  without  assuming  a  constant  I  ,  we 

out 

i  will  address  some  of  the  features  of  a  confocal  resonator 

which  are  significantly  different  from  those  of  the  Fabry- 


Reflectivities:  R^  R2 

Magnification:  M 

Gain  cavity  length:  t 

Gain  function:  g  (r,  z) 

Loss  function:  a 

Effective  reflectivity:  Reff  *  R^M"0 

Dimensionality:  n  *  1  (2D),  n  *  2  (3D) 

Common  focus  Z.  *  ML/(M  -  1) 

0  pi 

Threshold  condition:  Reff  exp  [  /Q  ( g  -  a)  dz]  s  1 

1  3l  + 

Rate  equations:  *  (g  -  a) 

r3&-  »„  -  *>  1ST  =  «*  +  “><zo  ‘  ■ 

Figure  1  Schematic  diagram  of  a  confocal  resonator  and  key  system 
equations. 


Figure  2.  Schematic  diagram  for  a  confocal  resi 
and  the  associated  Intensity  mesh-po 
(I)  and  oure  amplification  (II)  reoii 


Perot  resonators  and  the  factors  which  cause  the  complexity 
in  our  system,  as  follows: 

(i)  The  spatial  non-uniformity  of  the  gain  function. 

(il)  The  mirror  curvatures  wh£ch  guide  the  propaga¬ 
tion  of  the  forward  (I  )  and  backward  (I~) 
intensity  differently. 

(iil)  The  r-dependence  of  the  threshold  condition. 

(iv)  The  z-dependence  of+  the  total  intensity  (I++I~) 
and  their  product  I  I  . 

(r)  The  variation  of  the  output  intensity  in  the 
transverse  coordinate. 

We  note  that  in  the  Rigrod  analysis  for  a  Fabry-Perot  reso¬ 
nator  (with  M  -  1)  most  of  the  above  described  features  are 
ruled  out,  and  hence  simple  expression  for  the  efficiency  is 
available.  Our  numerical  results  of  the  rate  equation  will 
show  features  (iv)  and  (v).  It  is  seen  that  the  spatial 

variations  of  the  intensity  I-(z)  and  I  „(r),  which  are 

ou  z 

the  essential  feature  for  the  confocal  resonators,  cannot  be 
ignored,  e.g..  Equation  [1]  should  be  integrated  rather  than 
being  multiplied  by  a  constant  outcoupling  fraction  (1-M~n). 


III. 


ANALYTIC  TREATMENT 


For  tractable  results,  we  shall  first  investigate 
an  analytical  expression  for  the  resonator  efficiency  n. 
Assuming  a  constant  mean-intensity  in  the  gain  function,  the 
on-axis  threshold  condition  gives  a  simple  expression  for 
the  resonator  efficiency  (see  Appendix  A) 


no 

X 


X  r  80L 
g0L  LoqL  ~  X 

1/2  m  R;Jf, 
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C2] 

C33 


where  we  have  introduced  an  effective  reflectivity  R  -  - 

—  n  ®  1  a 

RgM  ,  R2  being  the  reflectivity  of  the  output  mirror  (for 
Rl  •  1 )  and  M  the  magnification  of  the  confocal  resonator 
with  dimension  parameter  n  -  1  for  2D  (cylindrical)  and  n  - 
2  for  3D  (spherical)  mirrors.  The  important  feature  of  the 
above  generalized  expression  is  twofold: 


e  For  a  Fabry-Perot  system,  M  -  1  Equation  [2] 

reduces  to  the  Rigrod  results;  and 

e  For  a  confocal  resonator  with  M  >  1,  Equation 
[2]  provides  a  new  expression  which  fits  the 
numerical  results  better  than  those  of  conven¬ 
tional  expression  with  a  coupling  fraction  (1  - 
M  )  . 

The  maximum  efficiency  (n*)  may  be  calculated  by 
the  optimizing  condition  dn/dX  *  0, 


with  the  associated  optimal  value 


80L  (/D  ~  D5 ' 


where  0  -  (g0L/oQL)  is  the  unsaturated  gain  and  loss  ratio. 
Some  of  the  salient  features  of  the  above  analytic  expres¬ 
sions  are: 


(i)  As  D  increases,  the  maximum  efficiency  increases 
and  reaches  unit  when  D  approaches  infinite; 

(ii)  In  a  given  reflectivity,  say  R.  -  R_  -  1 ,  the 
optimal  magnification  is,  from  Equation  [5], 
given  by 

"opt  ■  «P  c«oL  </5  -  5>3-  C6] 

which  is  shifted  to  the  larger  value  when  the  total  gain  of 
the  system  Increases  for  a  given  internal  loss.  He  should 
note  that  the  validity  of  Equation  [4]  is  limited  by  the 
inequality^ 


ft  -  (1  -  <  Cl  -  exp  (-  a0L)]/a0L  -  fg,  [7] 

Some  of  the  examples  are:  for  D  -  10,  f  •  0.467,  we  have 
f2  -  0.432  when  oQL  -  2;  for  D  -  20,  f1  -  0.6,  we  get  fg  - 
0.632  when  o.L  •  1.  Therefore  the  maximum  energy  extraction 


of  an  unstable  resonator  *  is  limited  by  the  cavity  length 
which  is  in  turn  limited  by  the  inequality,  Equation  [7]. 
For  systems  with  long  active  media,  the  maximum  efficiency 
is  given  by  the  smaller  quantity,  either  f1  or  fg.  We  shall 
compare  the  aforementioned  features  with  our  numerical  re¬ 
sults  later. 


IV.  NUMERICAL  RESULTS 

The  rate  equations,  (F.2)  in  Figure  1,  subjected  to 
the  steady-state  (on-axis)  condition  will  be  solved  numeri¬ 
cally.  Due  to  the  curvatures  of  the  mirror,  the  forward  and 
backward  fields  are  propagating  differently.  We  shall  in¬ 
vestigate  the  intensity  profiles  across  the  output  mirror 
and  the  propagation  of  the  forward  and  backward  field.  It 
will  be  shown  that  the  transverse  variation  of  the  Intensity 
across  the  output  mirror  is  mainly  caused  by  the  pure  ampli¬ 
cation  of  the  forward  intensity  when  it  enters  regime  (II), 
as  shown  in  Figure  2.  In  the  following  discussion,  we  shall 
consider  the  situation  where  |dl_/dr | << | dl_/dz |  and  assume  a 
constant  unsaturated  gain,  i.e.,  gQ  -  constant.  However  the 
total  gain  will  be  strongly  z-dependent  due  to  the  z-depen- 
dent  total  Intensity. 

In  order  to  employ  Newton's  method,  we  require  an 
explicit  form  of  the  total  intensity,  i.e.,  the  z-dependence 
of  the  gain  function  must  be  known.  As  shown  in  Appendix  B, 
in  contrast  to  the  planar-mirror  system,  I+(z)  and  I_(z)  and 
their  product  are  strongly  dependent  on  z  due  to  the  expan¬ 
sion  of  the  backward  field,  I_(z),  and  the  accompanying  gain 
saturation  which  in  turn  affects  the  z-dependence  of  the 
forward  Intensity,  I+(z).  The  numerical  procedures  are 
shown  in  Figure  3  with  greater  detail  shown  in  Appendix  B. 
The  Initial  intensities,  I+(z  -  0)  are  solved  from  the 
steady-state  condition  and  the  propagation  of  them  gives  us 
the  on-axis  intensity  profiles.  The  numerical  results  are 
shown  in  Figure  4. a.  The  off-axis  intensity  of  the  forward 
field  across  the  output  mirror  plane,  z  -  L,  is  shown  in 
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Figure  4  (A).  The  z-dependence  of  the  forward  (IT)  and  backward  (I”) 
intensity  Inside  the  confocal  resonator  with  magnifica¬ 
tion  M  »  (1)  2,  (2)  4  and  (3)  6.  System  parameters  are 
g0L  *  6  .  ,  a  ■  2  x  10"3  cm"1,  L  *  35  cm,  a  *  1.5  cm 

and  Ri  »  R,  *  1 . 


k.w  /;•  k.-;*cr  vjr-c w  jc'jc  nr.  •.■"jsr  bicv 


Figure  4 .  b .  it  Is  seen  that  stronger  r-dependence  is  asso- 
olated  with  smaller  magnif ication  and  the  intensity  profiles 
are  almost  linearly  dependent  on  the  transverse  coordinate, 
r.  This  feature  allows  us  to  evaluate  the  output  power. 
Equation  Cl],  in  a  simpler  way.  For  spherical  (2D)  mirrors, 
we  obtain 


»  . 
lv 

& 


P  -  wa2(1  -  m’2)  C1  ♦  »a3(1  -  M*3)  Cg/3, 


C8] 


for  an  output  intensity  given  by  IQut  -  C1  ♦  C2r.  It  is 

interesting  to  see  that  the  above  expression  reduces  to  the 

conventional  form  with  a  constant  outeoupling  fraction  (1  - 
•  2 

H  )  when  B  •  0.  We,  however,  should  also  note  that  both  A 
and  B  are  M-dependent  and  the  actual  explicit  form  of  P  in 
terms  of  M  are  not  analytically  available. 


II 
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DISCUSSION  AND  CONCLUSION 


He  summarize  the  numerical  results  In  the  following 
discussions : 


(1)  As  shown  in  Figure  5,  the  efficiency  based  on 
the  mean-field  approximation  (MFA),  Equation 
[2],  over-estimates  that  of  the  numerical  (ex¬ 
act)  results,  particularly  in  the  large  M  re¬ 
gime.  Moreover,  the  maximum  efficiency  is  also 
over-estimated  by  the  expression  of  Equation 
[4].  We  note  that  the  optimal  magnification 
based  on  the  near-field  information  does  not 
simultaneously  meet  the  requirement  of  good  beam 
quality,  i.e.,  high  focusability .  To  show  this, 
we  also  plot  the  M-dependence  of  the  center 
Intensity  of  the  far-field  in  Figure  5.  He  see 
that  high  efficiency  at  lower  M  and  high  beam 
quality  at  higher  M  are  two  competing  factors 
which  should  be  optimized. 

(ii)  The  effects  of  the  gain  and  loss  ration,  D  - 
gpL/a.L ,  on  the  efficiency  profiles  are  shown  in 
Figure  6.  It  is  easy  to  see  that  higher  D  val¬ 
ues  result  in  higher  efficiency  and  the  optimal 
magnification  is  shifted  to  lower  values  when  D 
increases.  These  features  based  on  the  numeri¬ 
cal  results  may  be  clearly  seen  from  the  approx¬ 
imate  expressions  Equations  [4]  and  [5]. 

(ill)  In  Figure  7,  we  show  the  effects  of  the  unsatu¬ 
rated  gain  for  fixed  values  of  D.  It  is  seen 
that  the  numerically  calculated  peak  efficien¬ 
cies,  depending  on  the  gain,  deviate  slightly 
from  those  of  the  approximate  expression,  Equa¬ 
tion  C 4 3 ,  in  which  nf  should  be  independent  of 
the  gain  as  far  as  the  D  value  is  fixed.  More¬ 
over,  the  general  feature  of  the  optimal  magni¬ 
fication  which  Increases  as  the  gain  increases 
is  well  described  by  Equation  [5]. 
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(iv)  In  Figure  8,  we  show  the  center  intensity  of  the 
far-field  which,  in  contrast  to  the  (near-field) 
efficiency,  is  a  strongly  increasing  function  of 
the  gain,  for  a  fixed  D  value. 

(v)  To  provide  some  guidance  of  the  optimal  design, 
we  plot  the  gain  and  magnification  diagram  in 
Figure  9.  For  a  given  magnification,  there  are 
two  values  of  the  unsaturated  gain  which  result 
in  the  same  efficiency,  say  55  percent.  Howev¬ 
er,  the  lower  one,  g. ,  is  preferred  as  far  as 
pumping  is  concerned.  We  note  that  the  maximum 
efficiency,  n*>  may  not  always  be  achieved  and 
is  restricted  by  the  values  of  M,  which  is  fi¬ 
nite.  For  example,  when  gQL  •  10,  efficiency  of 
55  percent  is  the  most  one  may  achieve.  To 
achieve  a  higher  efficiency,  one  may  "tune"  the 
gain  to  meet  the  optimal  condition  for  a  fixed 
magnification  or  vice  versa. 

(vi)  For  the  situation  that  some  of  the  system  param¬ 
eters  are  not  well-defined  or  may  not  be  measur¬ 
ed  accurately,  the  analytic  expressions  of  the 
output  power  and  the  related  optimal  condition 
are  very  useful  for  the  analysis  of  unknown 
parameters.  For  example,  the  measurements  of 
the  output  power  at  the  low  M  and  high  M  regimes 
together  with  the  beam  size  will  provide  us  the 
information  of  the  maximum  output  power  and  the 
corresponding  optimal  magnification.  For  great¬ 
er  detail,  see  Appendix  D. 


Figure  8.  Effects  of  the  unsaturated  gain  on  the  center  intensity  of  the  fi 
an  averaged  intensity  illuminates  on  the  output  spherical  mirror 
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APPENDIX  A 

Efficiency  In  the  Mean-Field  Approximation  (MFA) 
A  Generalized  Rigrod  Calculation 
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As  mentioned  previously,  most  of  the  Rigrod-type  calculations 
for  the  energy  extraction  efficiency  are  based  on  the  constant 
intensity  product  I+*I_  *  constant,  which  is  true  only  for  planar 
resonators.  For  unstable  resonators  with  curved  mirrors,  both 
1*1  and  (I  ♦  I  )  are  strongly  z-dependent  as  shown  in  our 
numerical  results.  In  the  following  discussion  we  shall  solve  the 
rate  equation  for  a  telescopic  resonator,  positive  branch,  confocal 
unstable  resonator,  within  a  mean-field  approximation(MFA) .  Analytic 
expression  for  the  energy  extraction  efficiency  will  be  derived  and 
the  limiting  case  with  Mai  (planar  resonator)  will  be  compared  with 
that  of  the  Rigrod's  results. 


For  a  resonator  with  uniform  gain  distribution  in  the  transverse 
direction,  the  2-0  rate  equation  is  reduced  to  the  following  1-0  case 


-kij— • 


(A .  1) 
(A. 2) 
(A. 3) 


which  has  the  formal  solution 


,  .  >2 


(A. 4) 
(A. 5) 


where  the  pre-exponential  z-dependence  factor  is  due  to  the  beam 
expansion  of  the  backward  field  in  a  2-0  spherical-mirror  system. 

We  readily  to  see  that,  from  (A. 4),  (A. 5),  the  intensity  product 


(A. 6) 


is  z-dependent  for  M>1,  and  is  a  constant  for  the  planar  case,  Mai. 

The  associated  gain  function,  for  a  homogeneous  broadening,  is  given  by 


_ & _ 

i  +vi- 


u 

1  cVit 


(A. 7) 
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Me  shall  now  solve  the  rate  equation 


There  are  several  ways  to  solve  the  rate  equation,  however, 
in  order  to  employ  the  MFA  for  the  total  intensity,,  we  substitute 
(A. 7)  into  ( A . 1 )  and  rearrange  it  in  an  integral  form 

[  1  ^  +  V* x-)]  (  a  .  9 ) 

■^1  6 

In  general,  the  above  equation  is  not  in  a  s  parated  form  of  z  and  1^ 
due  to  the  z-dependent  C2  and  the  Independent  total  intensity.  For 
analytic  form  of  these  integrals,  we  shall  impose  the  MFA  for  the 
total  intensity  in  the  RHS  of  (A. 9),  and  use  the  zeroth-order  of 
the  expansion  of  Cz  in  the  LHS  of  (A. 9).  Since  z/ZQ  <  1,  we  may  expand 

Cz  J.  (A. 10) 

Within  the  MFA, 

I+*)sJ, 

where  we  have  defined  the  mean  value 

and  we  have  used  the  boundary  conditions  in  (A. 8). 


(A. 11) 
(A. 12) 


Using  (A. 11)  and  (A. 12),  we  may  calculate  the  total  intensity 
in  the  MFA : 

-§ CA.13, 

catirA  €= *  /-AT  * 

which  may  be  re-written  as,  by  using  the  steady-state(on  axis)  condition 

e>  _  z 


RfQM  -  e 


(A. 14) 


+2€  l5) 


Within  the  MFA,RHS  of  (A. 9)  is  given  by 


(A. 16) 


We  shall  now  evaluate  LHS  of  (A. 9).  Using  the  leading  term  of 
(A. 10),  we  are  able  to  work  out  this  integral  follows: 

LHS 

(A'17> 


From 


LHS=RHS,  we  solve  for  the  intensity  on  the  output  mirror, I2  , 


UfAcrc. 


(A. 19) 


(A. 20) 


(A. 21) 


e  *  i- m~'  . 


WE  note  that  the  above  expression  reduces  to  that  of  Rigrod  for  the 
case  of  planar  mirrors  where  M*l,  and  hence  Xa-lntR^Rj),  with  <?*0. 
We  also  note  that  the  above  new  expression  for  the  output  intensity 
of  a  telescopic  resonator  shows  its  M-dependence  through  x,  f^,  f2 
and  €•  .  These  results  have  not  been  derived  in  the  literature, 
although  the  high  magnification  os* with  f^*l,  and  € *0,  has  been 
used  in  Ref.  2,3. 
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Knowing  the  output  intensity  acrossing  the  output  mirror,  I 2, 
we  are  able  to  evaluate  the  output  power  ,  from  Eq.(l), 

-wr-  ***o-*vz .  CA  22) 

where  we  have  used  a  mean  value  such  that  a  constant  outcoupling 
fraction  (l-M  )  is  defined  in  the  above  equation.  In  general  I2  is 
transverse  coordinate  dependent  and  (A. 22)  needs  numerical  method. 
Substituting  (A. 18)  into  (A. 22),  we  may  calculate  the  efficiency 
within  the  MFA  and  constant  outcoupling  approximation  as  follows. 

v "  ^ 'l ~jeF)  > 

/r-(^)f*^x).  <A-23) 


Note  that  2<F  is  the  first-order  correction  term  resulted  from 
the  expansion  of  Cz>  Eq.(A.iO).  In  the  zeroth-order  approximation 
and  when  R1aR2*1,»  EQ*(fl*23)  reduces  to  Eq.(2).  Furthermore,  as 
shown  by  our  numerical  results  that  the  output  intensity  is  only 
slightly  depending  on  the  transverse  coordinate  for  large  M  and 
a  constant  outcoupling  (l-M’2)  is  valid  in  Eq.(A.22).  The  overestimat¬ 
ion  of  Eq.(3)  may  be  caused  by  its  neglecting  the  correction  term. 

we  finally  note  that  the  MFA  expression  for  the  ef f iciency,  (4.2$), (2), 
fits  well  with  the  numerical  results  even  for  systems  with  high 
gain  and  loss  ratio,  within  5  X  deviation.  This  is  in  contrast  to 
that  of  the  Rigrod's  results  where  the  MFA  expression  deviates  from 
the  numerical  results  significantly  when  the  gain/loss  is  large. 
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We  shall  now  solve  the  rate  equation,  (A.l)  and  (A. 2)  subjected 
to  the  steady-state  condition,  we  first  should  point  out  that  the 
direct  numerical  solving  the  rate  equation  by,  e.g.,  Runge-Kutta  method 
is  difficult  due  to:  (i)the  boundary  conditions  in  the  regime  I  and. II 
and  (ii)the  intensity  dependent  gain  function.  In  order  to  solve  the 
rate  equation  by  Newton's  method,  the  explicit  form  of  the  gain  func¬ 
tion  must  be  known.  For  this  purpose,  the  total  intensity  is  imposed 

"  W+X9>*ZJZ (a  n 

which  is  numerically  justified  and  shown  to  be  valid  for  M<5. 


Using  (B.l),  the  gain  function  is  given  by 


(B.2) 


which  is  substituted  into  the  steady-state  condition  and  gives 


(B.3) 


where  may  be  sovled  by  Newton's  method  from  the  above 

self-consistent  equation. 

Given  the  explicit  form  of  the  gain, we  are  able  to  work  out 
the  forward  and  backward  intensities,  from  (A. 4)  and  (A. 5), 


(B.4) 


(B .  5) 


Above  expressions  give  the  on-axis  intensity  profiles.  For  intensity 
off-axis,  the  propagation  of  the  forward  field  must  be  treated  more 
carefully  in  the  regime  where  1  is  absent. 


As  shown  in  the  above  figure,  in  the  pure  amplification  regime  (X), 
the  gain  is  given  by 

9&m  1  +  *  (B.6) 

and  solution  of  the  rate  equation,  (A.l),  becomes 


(B.7) 


where  zQ  is  r-dependent  and  is  given  by 


(B.8) 


(B.7)  may  be  easily  worked  out  and  a  self-consistent  equation  is 
obtained  *  “i 


j  (B .! 


Again,  by  employing  Newton's  method,  we  are  able  to  solve  the 
output  intensity  acrossing  the  output  mirror  plane  for  a/M<r^  a. 
The  input  intensity  in  (B.9)  is  given  by  (B.4)  with  z*zQ  which  is 
now  r-dependent. 
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As  Show  In  the  above  figure,  the  far-field  amplitute  is  given  by 


where  is  the  first-order  Bessel  function.  Noting  that 
asf-+0,  the  center  intensity  of  the  far-field  is  given  by 

J^lUff-O))  =  2)  ,  (C . 2) 


ii- » 
*  ■. 


where  I_is  the  value  for  M*eo  and  we 
intensity  incident  upon  the  aperture 


tX~M  ] . 


have  used  the  mean  output 


( C .  3 ) 


(OKS 


we  note  that  I  t  depends  on  m 
(1-H*2)  is  an  increasing  function  o 
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and  the  outcoupling  fraction 
f  M,  therfore,  their  produt; 


I has  a  maximun  value  at  the  optimal  magnification  which  is  shifted 
to  larger  value  compared  with  that  of  the  near-field  efficiency.  For 
the  purpose  of  maximun  extraction  energy  as  well  as  high  focusability, 
we  shall  re-define  an  efficiency  based  on  the  far-field  energy  rather 
than  the  near-field  output  power  as  discussed  previously. 

The  new  efficiency  based  on  the  far-field  intensity  is  defined  as 

5 

C  icfl^L  '  (c,4) 

where  E^  is  the  center-ring  energy  of  the  far-field  given  by 

and  I(^)  is  the  far-field  intensity  given  by  |Up  in  Eq.  (C.l). 

We  note  that  in  an  unloaded  resonator  E^  increases  as  M  increases, 
however,  for  a  loaded  resonator  the  energy  extraction  is  saturated 
by  the  gain/loss  and  competing  with  the  outcoupling  fration  (1-M  ). 

Therefore, to  perform  a  resonator  with  high  energy  extraction  effici¬ 
ency  as  well  as  high  focusability,  we  shall  choose  an  optimal  magnif¬ 
ication  which  is  larger  far  larger  gain  systems. 


APPENDIX  D 

Analysis  of  the  Unknown  Parameters  via  the 
Measurements  of  Output  Power  and  Beam  Area 
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Due  to  the  high  cost  of  the  mirror,  particularly 
for  large  aperture  systems,  it  is  highly  unlikely  that  we 
may  change  the  system  magnification  continuously  to  obtain 
the  optimal  condition.  Therefore  it  is  highly  desirable  to 
estimate  the  optimal  magnification  by  some  of  the  measurable 
laser  parameters.  In  the  following  discussion,  we  shall 
suggest  some  possible  procedures  of  determining  the  unknown 
parameters  by  some  of  the  measurements. 

From  Equation  C 3 3 »  we  readily  see  that  experimental 
measurements  of  the  gain  and  loss  shall  enable  us  to  evalu¬ 
ate  the  optimal  magnification.  However,  when  the  measure¬ 
ments  of  the  gain  and  loss  are  inconvenient  or  inaccurate, 
there  is  an  alternative  to  evaluate  them.  To  demonstrate 
this  mathematically,  we  combine  the  expressions  of  Equations 
C 3 3  and  [4]  to  find  the  output  power 


out 
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As  shown  in  Figure  D,  the  output  power  has  slopes  m1  and  m2 
at  the  low  X  and  high  X,  respectively,  given  by 


■i  -  AIsX0/(o0L)' 


CD. 2a] 


m 


2 


[D.2b] 


where  XQ  -  (gQ  -  aQ)L  is  the  value  at  the  oscillation  thres¬ 
hold,  A  and  I  are  the  beam  area  and  the  saturation  intensi- 

9 

ty  at  which  g  -  gQ/2.  For  Equation  [D.2b],  we  see  that 
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x  [iklrf/Rfo)] 


Figure  D.  Schematic  diagram  of  the  output  power  versus  resonator  parameter 
X.  Note  that  measurements  of  the  slopes  and  will  provide 
us  the  optimal  value  X*.  * 


measurement  of  A  together  with  the  slope  m 2,  defined  by  the 
threshold  value  Xq  and  the  measurement  of  the  output  power 
closes  to  that,  enables  us  to  determine  I  .  From  the  mea- 

9 

sured  slope  m1 ,  XQ ,  A  and  the  calculated  I3>  we  then  may 
calculate  the  loss  which,  from  the  relation  XQ  ■  (gQ  -  oQ)L, 
gives  the  gain.  Knowing  the  gain  and  loss,  we  finally  may 
evaluate  the  maximum  efficiency  and  the  corresponding  opti¬ 
mal  magnification  from  Equations  [3]  and  [4]. 
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C ********  ?IITD  THE  OPTIMAL  KAGNIPICATIOK  POR  UHSTABLE  RESONATOR  ******* 

C *******  EXTRACT IOH  EPPICISHCT.  /RE. POR/  J.T.  LIE  2-10-84  ******* 

C ********  IS  THIS  PROGRAM  GAIH-GL(Z.O) ,  WITHOUT  R-DEPSHDBHCB  ********** 

C ********  RATE  EQUATION  IS  1-D  ,  BUT  ALLOW  IHTEHSITT  I.I(Y,R)  ********** 

1 ;  .ccttbr^par-pield  ( 2 ) 

IMPLICIT  RBAL*4(A-H,0-W,Z) 

DIMESSIOH  TT(2) ,XM1 (201 ) ,XM2(201 ) ,TC(201 ) ,TE(201 ) ,TR(201 ) ,IP(201 ) 
LOGICAL  ERROR 
BITERS AL  PT.GT 

COMMOI  2H,GL,AL,XM,R1  ,R2,XL,A,R,TT 

TTPS*, 'EETBR  GL,AL,XR,DH,HH«  5,0.25,1 . ,2.95,100* 

READ (5, 55)  GL, AL,XR,DH,HS 
55  P0RMAI(4G,I) 

C  INPUT  PARAMETERS:  R1 , 2(REPLECT ) ,  XL(CA7ITT  LESOTH),  A (OUTPUT  MIRROR  RADIU 
C  R( POSITIOS  IS  VERTICAL) ,  XM(MAGSIP) ,  (GL,AL)«( GAIN, LOSS) 

C  DIMESSIOH  0?  THE  RSSOIAIOR  SITES  BT  XS  XH«1(2D)  XH-2(3D) 

C .  refelctivity  of  the  output  lenses  at  193  am  wav*  (RX) 

RX-0.04 
R1-XR 
R2-1  . 

e . XL-cavity  length  only  affaeta  M*  But  not  efficiency,  her*, 

c . v«  w*  r*duo*  XL  By  a  factor  of  100.  to  avoid  underflow.... 

XL«1 . 

A*1 .5 
XS*2. 

DM  *DV /PLOAT ( IS ) 

C ******* ajoVE  hi  SITES  THE  HIGHEST  RAHGE  OP  XM,  PROM  1  TO  (RH+1 H3.95ETC) 
HPTS-HH+1 
DO  40  t«1 ,HPTS 
XM-1 .05  ♦  PLOAT (X-1 )*DM 
DO  33  1*1,2 

R«A*( FLOAT (1-1 )*{XM— 1 . )+1 . )/XM 

c.....  initial  value*  for  H*wton  method  z*1 .  ft  5*  for  PT  ft  GT 
X-1.0 

CALL  SEWT0I(X,7T, ERROR) 

C . above  result  of  X  will  Be  used  a*  the  input  intensity  in  GT 

TT-X 

1-5. 

CALL  HEWTOI(X,GT, ERROR) 

TT(I)-X 
33  COHTISUE 

CM1  *1 .  -  1  ./XM 
CM2-1 .  -  1 . /XH**2 
CM3-1.-1 • /XM**3 

c . ouput  intensity  approx.  By  linear  fuct  Betwwen  yt(1 )  ft  yt(2) 

C . on-axis  intensity  (TIC)-TT(I) . 

TIC-TTM  ) 

CB-(TT(2)-TT(1 ))/(CM1*A) 

CA*TT(2 )  -  CB*A 

C  THE  RBSOHATOR  EPPICIESCES  WILL  BE  HORMALIZSD  BT  GL*A**XS  POR  1ft2D(XH*1,2) 


RIGROD ■ ( GL/ ( AL  +0. 5*XH*AL0G(XM/SQRT(IR) ) )  -  1.) 

RD1  *(GL/(AL  +0. 5*ALC0(XM) )  -  1.) 

c .  to  oonpar*  exact(TE),  Rigrod(TR),  ooated(TC),  far  field  int(T?) 

TE(E)*(CA*CM2  ♦  ( 2 . *CB*A/3 • ) *CM3 ) /GL 
SR  ( X ) *0 . 5 *XH*ALOO ( XM/SQRT (XR ) ) *RIQR0D /GL 
C .  find  far-field  center  intensity,  with  seas  ist.(SA) . 


XA-CA+0.5*CB*A*(1 .+1 ./XM) 

T?(K)«EA*CM2**2/GL 
XN1 (I)-XM 
AO  COHTIHUB 

CALL  PLOTPT (XM1  ,YX,H?TS) 

CALL  PLOTPT  XM1  ,TR,HPTS 
CALL  PLOTPT (XK1, IF, HPTS) 
shp 

C *♦•***♦*  EXTERIAL  PUICTIOIS  ****** . 

PUHCTIOI  FT(X,FX,PPX) 

COMHO I  XH,GL,AL,XH,R1 ,R2,XL,A,R,TT 
C1-2.*X  +  1. 

C2-2.*X*XM  +  1. 

F1«AL0G(XM*C1/C2) 

F2-(  {UI-U  )/(2.*XM*OL)  )*(0L-AL-0. 5*ALOO(XH*n»/(Him2) ) ) 
FX-X  -  F2/F1 

PFX-1.  -  <P2/F1**2)*(XX-1.)/(C1*C2) 

RXTURH 

sip 

C 

FUXCTIOH  OT(I,FXf DPI) 

COMHOH  XH,GL,AL,XM,R1 ,R2,XL,A,R,TI 


C . find  intensity  I  (forward)  ia  ragiaw  where  I(back)-0 . , 

Z-XM*XL/ (XH>1 . ) 

XL«Z*(R/A  -  1 ./XX) 

ZZ-XL— XL 

T1«(1  —ZZ/(Z*(2.*TT+t . ) ) )**(2.*TT*Z*GL/IL) 

T-SQRI (R2)*tt*IXF( (GL-AL)*ZZ/XL)*T1 
01 «(  0L-AL-AI*X ) / ( OL-AIr-AL*T ) 

FX«AL00(AS3(X/T))  -  (OL>AL)*(XL/XL)  -  (OL/AL)*ALOO(ABB(01 )) 
PFX-1. /X  +  0L/(01*(0L-AL-AL*T)) 

RICTRf 

am  _  _  _ 

3UBRO0TIHB  HBMTOI  ( X ,  PUHC ,  XRROR  ) 

COMMOH  XH,0L,AL,XX,R1fR2,XL,A,R,Tr 
LOGICAL  ERROR 
ERROR-. FALSE. 

TOL-1 .OB-A 
MAX-50 

PO  20  1-1, MAX 
XI -X 

CALL  FUIC(X,FX,PFX) 

IP(PFX  .SQ.  0.0)  GOTO  99 

PX-FX/PFX 
X-X1  -  PX 

I?( ABS(DI)  .LE.  ABS(X*TOL) )  GOTO  30 
20  COHTIJUS 

XRROR— . TRUE • 

30  RSTURH 
99  ERROR-. TRUE. 

RXTURX 

ZIP 


END 
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